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Muon spin relaxation (/iSR) measurements in iron oxy-pnictide systems have revealed: 
(1) commensurate long-range order in undoped LaOFeAs; (2) Bessel function line shape in 
La(Oo.97Fo.o3)FeAs which indicates possible incommensurate or stripe magnetism; (3) anomalous 
weak magnetism existing in superconducting LaOFeP, Ce(Oo.84Fo.i6)FeAs, and Nd(Oo.88Fo.i2)FeAs 
but absent in superconducting La(Oo.92Fo.o8)FeAs; and (4) scaling of superfluid density and Tc in 
the Ce, La, and Nd-FeAs superconductors following a nearly linear relationship found in cuprates. 

PACS numbers: 74.90.+n 74.25.Nf 75.25.+Z 76.75.+i 



A renewed interest on high-Tc superconductivity has 
been generated by recent discoveries of iron-oxypnictide 
superconductors LaOFeP {Tc 5K) [l[ and LaOFeAs 
(Te - 26K) i, followed by subsequent development 
of materials with higer TcS up to ^ 55K containing 
rare earth (RE) elements, such as Ce, Nd, Sm 0, [1| 
instead of La. Carriers are doped by (0,F) or (La,Sr) 
substitutions as well as high-pressure oxygen syn- 
thesis For overall understanding of mechanisms of 
high-Tj, superconductivity, it would be very intstructive 
to compare these new superconductors with cuprate 
systems. Muon spin relaxation (/iSR) studies have 
provided unique information on magnetic order 0, [lo| 



and superfluid density 



11, m 



m cuprate 



systems. In this letter, we report /iSR measurements on 
undoped and doped iron-oxypnictide systems, LaOFeP, 
LaOFcAs, La(Oo.97Fo.o3)FeAs, La(Oo.92Fo.o8)FeAs, 
Ce(Oo.84Fo.i6)FeAs and Nd(Oo.88Fo.i2)FeAs. Our 
results demonstrate several generic features common to 
the iron-oxypnictides and cuprate systems, including 
long-range commensurate antiferromagnetism of the 
undoped parenet compounds, incommensurate or stripe 
magnetism of lightly doped systems near the border of 
magnetic and superconducting phases, and scaling of 
the superfluid density and Tc following a nearly common 
linear relationship. We also report anomalous weak 
magnetism detected in some of the superconducting 
systems. 

Polycrystalline specimens of all the FeAs compounds 
were synthesized at Beijing National Laboratory for Con- 
densed Matter Physics, following the method described 
in 15|. Specimen of LaOFeP was synthesized at Os- 



aka University, in the silica tubes using LaP, Fe and 
Fe203 as starting materials. These specimens were made 
into ceramic pellets of 8 mm in diameter and 1^2 
mm thick. The specimens of undoped LaOFeAs and su- 
perconducting La(Oo.92Fo.o8)FeAs are identical to those 
used in neutron scattering measurements [l6| which re- 
vealed collinear antiferromagnetic order of LaOFeAs be- 
low ~ 134 K and absence of long range magnetism in 
the superconducting La(Oo.92Fo.o8)FeAs. /iSR measure- 
ments were performed at TRIUMF using a He gas-flow 
cryostat in zero field (ZF), longitudinal field (LF), and 
weak-transverse-field (WTF) to characterize magnetism, 
and in transverse field (TF) to determine the magnetic 
field penetration depth. Details of the /kSR method can 
be found in refs. [TEIII- 

Figure 1(a) shows the time spectra of muon spin po- 
larization in zero field observed in the undoped parent 
system LaOFeAs. The long-lived and very clear muon 
spin precession signal indicates spatially long-ranged and 
homogenous magnetism, consistent with the commensu- 
rate Bragg peak for collinear antiferromagnetic structure 
found by neutron studies [l^, [l3| . The /xSR spectra can 
be fitted to a high frequency /(T ^ 0) ~ 23 MHz signal 
with a dominant (^ 60 % volume fraction) signal ampli- 
tude, and additional signal with much lower frequency 
/(T 0) - 3 MHz from a minor (~ 10% ) volume 
fraction. The temperature dependences of these frequen- 
cies are shown in Fig. 2(a). These results in LaOFeAs 
are essentially consistent with earlier reports by ^SR and 
Moessbaurer studies [l3,[l3| in LaOFeAs. Comparison of 
sub-lattice magnetization (cx f(T)) with spin-wave theo- 
ries will be published separately. 
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FIG. 1: (color) Time spectra of Zero-Field jiSK in (a) un- 
doped LaOFeAs and (b) 3%doped La(Oo.97Fo.o3)FeAs. The 
solid lines in (b) for T = 100 K represents a single Bessel 
function multiplied to a Gaussian damping, while for T — 2 
K a sum of three Bessel*Gaussian signals. 



ZF-/iSR spectra of the 3% doped La(Oo.93Fo.o7)FeAs, 
shown in Fig. 1(b), exhibit much faster damping than 
those in undoped LaOFeAs. The fast-damping spectra 
in Fig. 1(b) fit well to Bessel functions (multiplied by a 
Gaussian damping factor) , which are shown by the solid 
lines. The Bessel function line shape was first observed in 
ZF-^SR in an incommensurate spin density wave (SDW) 
system TMTSFzPFg [20[. Subsequently, this line-shape 
was also observed by ZF-/iSR in Lai.875Bao.i25Cu04 [2l|, 
Lai.47Ndo.4Sro.i3Cu04 and several other cuprate systems 
[lo| which undergo formation of static spin stripes. The 
comparison of parent and doped systems in Fig. 1 ex- 
hibits exactly the same features as in the case of the 
cuprates, La2Cu04 and the 1/8 doped spin stripe sys - 
tems, shown in Fig. 4(b), 2(a) and 2(b) of ref. [l0| . 
On the other hand, it is in principle possible to expect a 
highly damped spectra also in a commensurate antifer- 
romagnet with substantial randomness, such as electron- 
doped cuprate systems in the antiferromagnetic region 
[23 . [2^. Neutron scattering studies are required for a 
clear distinction of these three spin structures and deter- 
mination of the spatial periodicity. 

Figure 2(a) compares the ZF precession frequencies 
of the two cosine signals in undoped LaOFeAs (closed 
and open red symbols) with the frequencies of the Bessel 
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FIG. 2: (color) (a) Muon spin precession frequency ob- 
served in zero field in LaOFeAs (for two cosine signals) and 
La(Oo.97Fo.o3)FeAs (for three Bessel signals), (b) The mag- 
netic volume fraction estimated by WTF measurements with 
the external field of 50 G. 



function signals in the 3% -doped compound (closed blue 
symbol). In the case of incommensurate SDW order, the 
frequency of the Bessel function corresponds to the in- 
ternal field at the muon site near the maximum SDW 
amplitude. The lower Bessel frequencies at T = 2 K in 
the 3%doped system implies that the average internal 
field at the muon site is substantially reduced from that 
in the undoped LaOFeAs. This is remiscent to the case of 
cuprates, where the ZF-/iSR frequency in La2Cu04 was 
about '--^ 30 % higher than the Bessel frequency in the 1/8 
doped stripe systems, as shown in Fig. 3(b) of ref. [lO|. 
Such a large variation of internal field via 3% carrier dop- 
ing is hardly expected by a simple magnetic dilution of 
commensurate antiferromagnets. Although /iSR is a lo- 
cal magnetic probe which cannot determine precise spin 
structures, these considerations suggest a strong possi- 
bility of incommensurate or static stripe magnetism in 
La(Oo.93Fo.o7)FeAs. 

Figure 2(b) shows the volume fraction of magnetically 
ordered region in the undoped and doped FeAs systems, 
obtained in the WTF-/iSR measurements [24,] with WTF 
~ 50 G, after a correction for a background signal with 
the spectral weight of less than 10% of the total signal 
amplitude. In both doped and undoped systems, the 
magnetic order develops in essentially full volume frac- 
tion below the Neel temperture Tjv. A continuous change 
with finite volume fraction was observed in a narrow 
range of temperature within 4-/- 5 K of Tjv in both sys- 
tems, indicating a weakly first order nature of thermal 
phase transition. A rather high Tjv 110 K and full 
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FIG. 3: (color) Muon spin relaxation function observed in zero 
field in (a) La and Ce-FeAs superconductors, (b) Nd-FeAs and 
LaOFeP superconductors, and (c) in the parmagnetic state of 
LaOFeAs and La(Oo.97Fo.o3)FeAs. 



magnetic volume fraction of the doped system indicate 
robustness of static magnetism against carrier doping, 
which is reminiscent to the case of non-superconducting 
electron-doped cuprates. 

ZF-/iSR measurements have also revealed an anoma- 
lous and weak relaxation existing in normal and su- 
perconducting states of LaOFeP, Ce(Oo.84Fo.i6)FeAs 
and Nd(Oo.88Fo.i2)FeAs, as well as in the paramag- 
netic state of LaOFeAs and La(Oo.97Fo.o3)FeAs above 
T/v- As shown by the time spectra in Fig. 3(a)- 
(c), this relaxation is seen in most of the systems 
in the present study, except for the superconducting 
La(Oo.92Fo.o8)FeAs. The anomalous relaxation has also 
been reported in SmOo.82Fo.i8FeAs [25[, while the ab- 



sence of the effect in the La-based superconducting sys- 
tem is consistent with an earlier report [1^. These ZF- 
/iSR spectra exhibit almost no temperature dependence, 
except for the Ce compound below T ~ 4 K and Nd 
compound at T = 2.4 K, which is presumably due to im- 
minent ordering of Ce or Nd moments. The static origin 
of this anomalous relaxation was confirmed by LF-/zSR 
measuremens in LaOFeP at T = 8 and 2K, while static 
and additional dynamic effects were found for the Sm 
compound in ref. [2^ at T = 60 K. The observed relax- 
ation rate of 0.3 - 1 /.ts"^ corresponds to about 3 - 10 G 
of random static internal field. If it comes from dilute 
frozen moments of ~ 1 Bohr magneton, this field corre- 
sponds to 0.01 - 0.1 % level of concentration per formula 
unit. These features indicate: (1) the slow anomalous 
relaxation does not necessarily correlate with supercon- 
ducting transition; (2) this effect exists not only in sys- 
tems containing magnetic rare earth elements but also in 
several La based systems; (3) dilute frozen moments, con- 
ceivable for a non-stoichiometry of Fe and/or rare-earth 
sites and/or for minority impurity phases, are possible 
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FIG. 4: (color) Muon spin relaxation rate a observed in 
transverse external fields in superconducting FeAs and FeP 
based systems. The relaxation envelope was fit to a pre- 
fixed exponential function multiplied by a Gaussian function 
ea;p(— (T^t^/2). Absence of field dependence is consistent with 
the effect due to the magnetic field penetration depth. 



origin of the observed effect, but further studies are re- 
quired for conclusive determination. 

On the four superconducting specimens (those in Figs. 
3(a)and (b)), we performed TF-/iSR measurements in TF 
= 300, 600, and 1200 G to measure the magnetic field 
penetration depth A, and to derive the superfluid den- 
sity ns/m* (superconducting carrier density / effective 
mass) proportional to the relaxation rate a due to flux 
vortex lattice below Tc- To correct for the effect of the 
anomalous magnetic relaxation shown in Fig. 3, we first 
obtained an exponential relaxation rate in TF above Tc, 
fixed this value in the fitting, multipled a Gaussian pre- 
cession envelope exp{—a'^t^/2) to this exponential damp- 
ing, and derived a as shown in Fig. 4. This procedure 
does not make essential difference from a simpler fit- 
ting with Gaussian damping alone in the FeAs based su- 
perconductors, since the pre- fixed exponential relaxation 
rate (0.03, 0.24 and 0.34 /xs~^, respectively, for the La, 
Nd, and Ce compound) is much smaller than the super- 
conducting Gaussian relaxation rate. For LaOFeP, the 
exponential (1.36 /is~^) and Gaussian relaxation rates 
are very close, resulting in a significant systematic uncer- 
tainty in accuracy of the superfluid density. The rapid 
increase of a in the Nd compound (at T=2K) and Ce 
compound (below T = 4 K) are presumably due to im- 
minent magnetic order of moments on these elements. 
Otherwise, a{T) shows a saturation at low temperatures, 
consistent with an earlier report [2^ . In view of difficulty 
in determining the pairing symmetry using /j,SR results 
on ceramic specimens experienced in the cuprates, how- 
ever, we postpone the symmetry arguments until results 
on single-crystal specimens become available. 

In the cr(T 0) vs. plot of Fig. 5, we compare 
the results of the FeAs based systems to cuprates and 
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ical work [30|. In this model, superconductivity should 
be absent in the doping region with static Cu/Fe mag- 
netism (except for phase separation), which is the case 
both in the copper oxide and FeAs superconductors. In- 
elastic neutron and two-magnon Raman studies of J and 
hujAF of FeAs systems would be quite interesting from 
this view point. In summary, we demonstrated several 
common features between cuprates and oxy-pnictides by 
comparing /j,SR results on their magnetism and super- 
fluid density. 
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973 project of MOST (China) at lOP, Beijing. 



FIG. 5: (color) A plot of the muon spin relaxation rate 
a{T —f 0) CK ns/m* versus Tc for the iron oxypnictide sys- 
tems (colored solid circles for the present results and colored 
open circles for ref. f26| compared with the results for the 
cuprates and alkali-doped Ceo systems ( 11, 28, 29]). The 
FeAs superconductors follow the linear relationship found for 
underdoped cuprates in this plot. 



a few other exotic superconductors. Three points from 
the present study and two other points from ref. (26j 
demonstrate that the electron-doped FeAs systems fol- 
low a nearly linear relationship between Tc and ng/m* cx 
(t(T — > 0), with the slope approximately same as that 
found for cuprates. We note that our earlier results on 
an electron-doped cuprate (Nd,Ce)2Cu04 
low the behavior of hole-doped cuprates. 
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also fol- 
Figure 5 de- 

omnstrates that many type-II superconductors, including 
cuprates, iron-oxypnictides, and Ceo systems, share ap- 
proximately same ratios between their superfluid energy 
scales and Tc- This implies that Us/m* is an important 
factor in determining Tc- 

A large relaxation rate cr(T — > 0) observed in a rather 
low Tc LaOFeP, however, suggests that additional fac- 
tor(s), such as closeness to the magnetically ordered 
state, also plays a significant role in determining Tc, as 
discussed in ref. J2», '29*1 for cuprates. In these refer- 
ences, one of us [28i, i29|J proposed a paring mechanism 
for cuprates based on a charge motion (with an effec- 
tive Fermi energy kTp derived from the superfluid den- 
sity) resonanting with antiferromagnetic spin fluctura- 
tions (with an energy scale Hloaf represented by the ex- 
change interaction J) as kTp Tiljaf- With the present 
results suggesting (t(T ^ 0) cx ns/m* cx kTp (in 2-d 
sysems) cx Tc, this resonant spin-charge motion model 
might also be applicable to oxy-pnictides if their Tiloaf is 
approximately 50 meV or so, as estimated in a theoret- 
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